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1.  Introduction. 

Logistics  research  covers  many  different  areas  related  to  problems  of 
designing  and  controlling  large  dynamical  systems.  In  recently  published 
books  [3,6]  one  can  find  several  good  discussions  of  these  problem  areas  and  the 
modern  trends  in  logisitcs  research.  This  research  is  carried  for  the  pur- 
pose of  developing  efficient  management  procedures  of  the  total  support 
functions  of  large  systems.  We  include  today  among  these  functions  not  only 
the  traditional  ones,  namely:  design  of  facilities;  production  planning  and 
control,  assignment  and  allocation  problems,  inventory  control,  transporta- 
tion and  others,  but  also  problems  areas  concerned  with  the  reliability  and 
maintainability  of  weapon  systems;  tactical  mini-  and  maxi-combat  problems, 
etc.  In  all  the  areas  mentioned  above  statistical  research  occupies  a central 
position.  One  can  find  hundreds  of  papers  that  have  been  written  on  various 
stochastic  and  deterministic  optimization  methods  in  each  one  of  the  areas 
mentioned  above.  Not  much  is  available,  however,  in  the  literature  on  sta- 
tistical theories  and  methods.  This  is  partly  due  to  the  complexity  of  the 
problems  and  partly  to  the  widespread  practice  of  applying  the  established 
"optimal"  procedures  in  which  simple  estimates  of  the  unknown  parameters  are 
substituted.  In  many  cases  such  a simple  approach  yields  adequate  results. 

It  cannot,  however,  be  always  recommended,  as  will  be  discussed  later. 

The  objective  of  the  present  paper  is  to  expose  and  discuss  some  sta- 
tistical problems,  which  the  author  has  studied  in  the  past,  and  to  indicate 
fruitful  problem  areas  for  further  research.  In  our  exposition  we  will  dis- 
cuss the  following  problems: 

1)  Statistical  control  of  two-echelon  multi-station  inventory 
systems; 

2)  Adaptive  manpower  forecasting  model; 

3)  Strategies  of  crossing  mine-fields. 


i 
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These  problems  areas  are  not  the  only  interesting  ones.  There  are 
other  interesting  statistical  problems  in  the  field  of  logistics  research, 
like  problems  of  evaluating  the  readiness  of  systems,  reliability  estima- 
tion, sequential  determination  of  surveillance  epochs  and  others.  For  some 
discussion  of  these  problems  see  Chapter  10  of  Modern  Trends  In  Logistics 
Research  [ 6]. 

2.  Two-Echelon  Multi-Station  Inventory  Problems. 

In  the  present  section  we  give  an  example  of  a complicated  inventory 
control  problem  which  requires  statistical  techniques,  since  the  information 
on  the  demand  distributions  is  incomplete.  The  problem  we  consider  is  that 
of  determining  optimal  stock  policies  in  a two-echelon  multi-station  model, 
designed  for  naval  applications.  The  two  echelons  under  consideration 
are  the  depot,  D (upper)  and  tender-ships,  T^,...,T^  (lower).  The  customers 
arriving  at  the  stations  are  submarines.  We  assume  that  the  monthly  demand 
for  a specified  item,  at  station  T^  is  a random  variable  , j = l,...,k. 


Fig.  1.  Inventory  flow  in  a two- 

echelon  multi-station  model. 
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The  flow  of  material  in  the  system  is  described  in  Fig.  1.  The  stations  can 

order  from  the  depot,  D,  at  the  beginning  of  the  n-th  month  a quantity  Y.  , 

J .n 

or  could  send  back  to  the  depot  unneeded  stock.  If  Q . denotes  the  stock 

n,  j 

level  at  T,  at  the  beginning  of  the  n-th  month  then  Y . > -Q  , , for 
J ^ n»J  n.j 

all  j * l,...,k.  Moreover,  2 Y+  . < V , where  Y+  “ max(0,Y)  and  V 

n,j  — n n 

designates  the  stock  level  at  the  depot  at  the  beginning  of  the  n-th  month. 

The  depot  orders  from  the  manufactures  at  the  beginning  of  the  n-th  month  a 

quantity  Z^,  which  arrives  at  the  depot  only  1 months  later.  The  lead- 

time  between  the  upper  and  the  lower  echelons  is  one  month.  The  objective 

is  to  determine  at  the  beginning  of  each  month  order  levels  (Z  ,Y  ,,..., Y ,) 

n n , I n, k 

so  that  the  expected  total  lower  echelon  cost 

k + 

C(Q  ,X  ) - 2 (c . (Q  - X ) + p (Q  - X ) } 

^ ^ n > i i n ,1  n f 1 

will  be  minimized  and  the  probability  that  at  each  month  the  depot  will  be 
able  to  satisfy  the  demand  of  the  stations  will  not  be  smaller  than  a pre- 
assigned tolerance  level  y.  Zacks  and  Fennell  [11,12]  studied  the  problem 
when  only  partial  information  on  the  distributions  of  the  demand  variables 

X .,..., X , is  available.  Bayes  adaptive  control  procedures  were  developed, 
n , i n , K 

These  procedures  determine  after  every  months,  according  to  the  total  observed 
demand  at  each  station,  the  predictive  distribution  of  the  future  demand 
variables.  The  predictive  distributions  obtained  from  the  Bayesian  model  are 
negative  binomial  with  adaptively  changing  parameters.  The  negative  binomial 
distributions  were  previously  fitted  by  Haber  and  Sitgreaves  [ 4 ] to  the 
demand  for  repair  parts  in  the  Polaris  fleet.  Such  distributions  were  also 
fitted  by  Denikoff,  Fennell,  Haber  and  Solomon  [2],  It  was  shown  in  these 
previous  studies  that  the  negative  binomial  distributions,  which  depend  on 
two  parameters,  can  provide  an  adequate  fit  to  the  demand  distributions  of 


most  classes  of  Items  in  the  submarine  fleet.  In  the  Bayesian  model,  if  we 


i 
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assume  that  demand  variables  X^X^...,  for  a certain  part,  at  different 

months  are,  conditionally  an  0 , independent  and  identically  distributed 

having  a Poisson  distribution  with  monthly  average  0 , and  if  0 has  a 

prior  ganrna  distribution,  G(—  ,v) , then  the  expected  distribution  of  X , 

T n+1 

n 

given  that  Z)  X = T , is  the  negative  binomial  distribution  with 
i-1  1 


p.  d. f . 


P1Vi ' ‘I**  • 1 


r(v  + i)  . v 

n J r(i+l)r(v  ) *n(1  " V * 1 0,1,, 


with  <|> 


, . , . , . and  v - v + T 
n 1 + (n+l)x  n n 


This  distribution  is  called  the  predictive  distribution  of  X given  the 

n+l 

past  history  of  the  demand.  The  model  is  thus  much  more  flexible,  since  the 
predictive  distributions  applied  for  the  determination  of  the  optimal  stock 
levels  are  not  the  same  every  month,  but  keep  adapting  to  the  actual  results. 
Algorithms  for  the  determination  of  the  optimal  order  levels  at  the  lower 
echelon  and  the  upper  echelon  were  developed.  For  the  proof  of  optimality 

of  the  overall  procedure  see  Zacks  [8],  The  cost  function  looks 

simple  enough,  but  due  to  the  structure  of  the  inventory  flow,  the  determi- 
nation of  the  order  vector  should  minimize  the  expectation 

k + + 

R(Y  ;Q  ) - S E {c  [ (Q  . + Y - X ) - X . .] 

n n i**l  i n § i n j 1 n j i n+ 1 > 1 


This  is  the  expected  cost  due  to  shortage  or  excessive  stock  at  the  end  of 

the  (n+l)st  month,  computed  according  to  the  predictive  distribution  of 

X . and  X . , .,  determined  at  the  end  of  the  (n-l)st  month.  The  program 
n , l n+l , l 

for  the  determination  of  the  optimal  Y is  not  too  complicated  and  the 
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computations  are  quite  fast.  On  the  other  hand,  the  determination  of  the 

optimal  upper  echelon  ordering,  Z , for  a lead-time  of  i • 2 months 

n 

requires  to  find,  at  the  end  of  the  (n-l)st  month  the  smallest  integer  z 
which  satisfies  the  inequality 


P{z  + Z , + V > 
n-1  n — 


k 

2 

i-1 


(Yn,i  + Vl,i> 


k 

+ S 

i-1 


(k 


n+2,i 


Qn+2,i)  ^n-l} 


where  (k^^  j “ ^n+2  i^  is  the  desirable  ordering  level  of  the  i-th  sta- 
tion two  months  in  the  future.  This  conditional  probability,  which  depends 
on  the  total  observed  demand,  at  all  the  k stations  during  the 

first  n-1  months,  is  a very  complicated  predictive  distribution.  In 
order  to  approximate  it  by  simpler  predictive  distributions  we  have  per- 
formed extensive  simulations  of  such  systems,  with  the  exact  computations 
of  the  Z^  values.  In  the  following  table  we  illustrate  a Six-Station 
case  where  the  demand  variables  at  each  station  are  independent  Poisson 
with  different  means,  X. 


Table  1.  The  parameters  of  a 6-station  system. 


i 

Vi 

Ti 

«i 

■ 

ci 

Pi 

Xi 

V 

Z 

1 

19. 

2. 

0 

2. 

12. 

2. 

- 

- 

2 

19. 

2. 

2 

4. 

10. 

4. 

- 

- 

3 

19. 

2. 

4 

6. 

8. 

6. 

- 

- 

4 

19. 

2. 

6 

6. 

8. 

8. 

- 

- 

5 

19. 

2. 

8 

4. 

10. 

10. 

- 

- 

6 

19. 

2. 

10 

..  . - 

2. 

12. 

12. 

- 

- 

Depot 

- 

_ 

- 

- 

1 

0 

wmar~ 
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Table  2.  A 12-months  simulation  of  a 6-station  system. 


n 

B 

□ 

B 

Y2 

Y3 

Y4 

Y5 

Y6 

V 

Z 

mm 

0 

2 

4 

6 

8 

10 

5 

0 

0 

0 

0 

0 

5 

62 

108 

kb 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

314 

Wm 

0 

0 

0 

0 

0 

0 

10 

10 

7 

10 

11 

14 

62 

63 

94 

Bfl 

9 

6 

2 

1 

0 

0 

0 

0 

0 

0 

0 

7 

7 

12 

312 

iM 

9 

0 

0 

0 

0 

0 

0 

10 

9 

10 

15 

19 

63 

66 

86 

mm 

6 

2 

0 

1 

3 

4 

0 

0 

0 

0 

0 

12 

12 

19 

252 

:K9 

4 

0 

0 

0 

0 

4 

3 

10 

10 

12 

15 

16 

66 

68 

120 

8 

5 

7 

2 

4 

7 

9 

0 

1 

0 

3 

6 

9 

19 

19 

98 

3 

A 

0 

0 

2 

14 

5 

8 

13 

15 

16 

8 

68 

65 

202 

7 

9 

11 

4 

10 

11 

0 

0 

0 

9 

5 

8 

22 

21 

252 

6 

4 

9 

5 

7 

6 

2 

7 

2 

10 

11 

17 

65 

49 

200 

5 

9 

0 

8 

8 

17 

2 

2 

12 

7 

9 

5 

37 

38 

156 

This  table  show  that  after  a few  months  of  correcting  a system  which  is  ill 

conditioned,  the  system  stabilized,  and  the  optimal  Z values  are  approxi- 

m 

mately  equal  to  the  sum  of  the  values  only  for  several  months.  We  have 

also  shown  in  [11]  that  the  Y and  Z values  which  are  determined  by  using 

n n jo 

Poisson  distributions  with  estimates  of  the  X values,  yield  average  monthly 
cost  larger  than  the  ones  obtained  by  the  Bayes  adaptive  procedure.  This 
illustrates  well  the  need  for  the  development  of  proper  statistical  procedures. 

3.  Adaptive  Manpower  Forecasting  Models. 


In  the  present  example  we  would  like  to  present  a forecasting  problem 
which  does  not  belong  to  the  class  of  problems  which  are  suitable  to  the 
commonly  used  time-series  forecasting  like  those  discussed  by  Box  and  Jenkins 
f 1 J . The  problem  is  that  of  periodic  forecasting  of  the  total  size  of  the 
Marine  Corps,  and  the  size  of  certain  parts  of  the  force.  The  need  for  a 
statistical  method  of  forecasting  the  size  of  the  force,  despite  the  fact 
that  individuals  belonging  to  the  force  sign  service  contracts,  is  due  largely 


to  the  problems  of  random  attrition.  This  attrition  is  not  of  constant 
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proportions  at  every  time  period  and  at  different  units  of  the  force.  With- 
drawal from  service  among  the  career  force  personnel  is  of  different  char- 
acteristics than  that  of  first  enlistees.  The  factors  which  significantly 
affect  the  attrition  of  the  first-term  enlistees  were  studied  by  Haber  [4], 

It  was  found  that  among  many  different  possible  factors  the  ones  which  sig- 
nificantly influence  the  retention  rates  are:  length  of  contract,  year  of 

entry  to  the  force,  education  and  race.  These  factors  are  not  independent; 
moreover,  the  year  of  entry  to  the  force  is  an  important  variable  because 
many  of  other  possible  factor  of  apparent  theoretical  importance  (economical, 
sociological,  psychological  and  political)  are  actually  time  dependent. 

Thus,  if  we  classify  the  individuals  according  to  the  year  and  quarter  of 
entry,  length  of  contract,  education  (less  than  high  school,  at  least  high 
school)  and  race  (white  and  non-white)  all  the  other  factors  do  not  contri- 
bute significantly  to  the  improvement  of  the  forecasting.  On  the  other 
hand,  different  parts  of  the  career  force  show  high  similarity  of  retention 
rates.  It  was  found  that  one  retention  rate  can  be  used  for  the  whole 
career  force  at  each  forecasting  period.  The  career  force  rate  shows  also 
small  variability  over  time. 

The  methodology  of  forecasting  applied  in  our  studies  is  based  on  the 
estimation  of  the  retention  probabilities  of  certain  subpopulations  over 
the  last  quarter,  and  determining  prediction  (or  tolerance)  intervals  having 
the  characteristic  that,  with  confidence  probability  (1-a)  we  anticipate 
the  size  of  the  groups  under  consideration  to  belong  at  the  end  of  the  next 
quarter  to  the  intervals.  Several  methods  of  determining  the  prediction 
intervals  were  compared. 


mm 
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In  the  following  table  we  illustrate  the  prediction  intervals  determined  by 
three  methods  (conditional  maximum  likelihood  - CMLE;  tolerance  limits  and 
Bayes  prediction  limits)  as  determined  on  actual  Marine  Corps  data. 


Table  3.  Prediction  limits  for  the  size  of  cohorts  on  July  2,  1972 
(a  = .05).  Forecasting  on  January  1,  1972. 


i 

j 

k 

m 

CMLE 

Tolerance 

Bayes 

Actual 

1 

1 

1 

1 

1401. 

1434. 

1382. 

1446. 

1393. 

1439. 

1413. 

1 

1 

2 

1 

507. 

520. 

496. 

523. 

502. 

521. 

499. 

1 

2 

1 

1 

3272. 

3303. 

3253. 

3315. 

3264. 

3307. 

3280. 

1 

2 

2 

1 

676. 

693. 

665. 

697. 

671. 

694. 

674. 

2 

1 

1 

1 

935. 

979. 

912. 

997. 

926. 

986. 

1009. 

2 

1 

1 

2 

876. 

916. 

855. 

932. 

867. 

922. 

910. 

2 

1 

1 

3 

1133. 

1165. 

1115. 

1177. 

1125. 

1170. 

1155. 

2 

1 

2 

1 

113. 

130. 

103. 

136. 

109. 

132. 

135. 

2 

1 

2 

2 

98. 

111. 

89. 

114. 

94. 

112. 

104. 

2 

1 

2 

3 

158. 

169. 

149. 

171. 

154. 

170. 

163. 

2 

2 

1 

1 

1973. 

2009. 

1952. 

2023. 

1964. 

2014. 

2019. 

2 

2 

1 

2 

2277. 

2309. 

2257. 

2322. 

2269. 

2314. 

2301. 

2 

2 

1 

3 

1112. 

1133. 

1099. 

1140. 

1106. 

1136. 

1121. 

2 

2 

2 

1 

246. 

260. 

235. 

264. 

241. 

261. 

256. 

2 

2 

2 

2 

223. 

236. 

214. 

239. 

219. 

237. 

230. 

2 

2 

2 

3 

171. 

181. 

163. 

183. 

167. 

182. 

177. 

3 

1 

1 

1 

1467. 

1518. 

1440. 

1540. 

1456. 

1527. 

1570. 

3 

1 

1 

2 

1692. 

1749. 

1663. 

1773. 

1630. 

1758. 

1781. 

3 

1 

1 

3 

2261. 

2325. 

2223. 

2352. 

2247. 

2335. 

2347. 

3 

1 

1 

4 

2554. 

2616. 

2522. 

2643. 

2541. 

2627. 

2633. 

3 

1 

1 

5 

2952. 

3001. 

2924. 

3022. 

2940. 

3010. 

2979. 

3 

1 

2 

1 

147. 

166. 

136. 

173. 

142. 

169. 

165. 

3 

1 

2 

2 

167. 

187. 

156. 

194. 

162. 

190. 

190. 

3 

1 

2 

3 

266. 

291. 

252. 

300. 

260. 

294. 

292. 

3 

1 

2 

4 

279. 

300. 

266. 

207. 

274. 

302. 

294. 

3 

1 

2 

5 

391. 

404. 

380. 

403. 

386. 

405. 

386. 

3 

2 

1 

1 

3333. 

3379. 

3307. 

3399. 

3322. 

3387. 

3403. 

3 

2 

1 

2 

4623. 

4676. 

4594. 

4699. 

4611. 

4685. 

4675. 

3 

2 

1 

3 

3380. 

3427. 

3354. 

3447. 

3369. 

3435. 

3416. 

3 

2 

1 

4 

4115. 

4159. 

4090. 

4177. 

4104. 

4166. 

4108. 

3 

2 

1 

5 

3207. 

3240. 

3187. 

3252. 

3199. 

3245. 

3207. 

3 

2 

2 

1 

331. 

348. 

319. 

353. 

325. 

350. 

345. 

3 

2 

2 

2 

431. 

452. 

417. 

459. 

425. 

455. 

451. 

3 

2 

2 

3 

411. 

432. 

397. 

440. 

405. 

435. 

428. 

3 

2 

2 

4 

466. 

488. 

453. 

496. 

460. 

491. 

478. 

3 

2 

2 

5 

526. 

540. 

516. 

544. 

522. 

541. 

526. 

Sums 


48170  49249 


47537  49656 


47902  49404 


49130 
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In  the  above  table,  the  index  i (1=1,2, 3)  designates  the  three  types  of 
contracts  (2,3,4  yrs).  The  index  j (j=l,2)  designates  the  race  (W,N-W) 
The  index  k (k=l,2)  designates  the  education  (LHS,AHS)  and  the  index 
m designates  the  time  of  entering  to  the  force  according  to  the  following 
code: 


i 

m 

Entry 

1 

i 

Jan.  1971 

2 

i 

Jan.  1970 

2 

July  1970 

3 

Jan.  1971 

3 

1 

Jan.  1969 

2 

July  1969 

3 

Jan.  1970 

4 

July  1970 

5 

Jan.  1971 

If  we  denote  by  t*le  size  of  the  (i,j,k,m)-th  subpopulation 

after  t periods  (of  six  months)  then  the  basic  theoretical  model  is  that, 

for  each  such  subpopulation  the  conditional  distribution  of  X...  (t)  given 

ijkm 

Xijkm(t_1)  is  the  binomial,  B(XiJkn(t-l) , 9iJkin(t));  where  eiJtan(t>  is 
the  retention  probability  for  that  subpopulation  during  the  t-th  period. 


The  CMLE  method  provides  (l-a)-prediction  intervals  for  X,.,  (t),  given 

ij  km 

Xijkm^t-^’  w^en  retention  probabilities  are  estimated  from  the  data 
by  the  conditional  maximum  likelihood  estimators.  The  tolerance  limits 
and  the  Bayes  limits  are  constructed  by  considering  first  the  2ARCSIN(/9) 
transformation; 


Yijkm(t)  “ 2810 


-1 


'WILL5 


This  transformation  has  desirable  characteristics  in  terms  cf  variance  sta- 
bilization even  for  medium  size  samples.  We  will  not  provide  here  the  de- 
rivations and  the  formulae  of  the  prediction  limits.  The  interested  reader 


-10- 


1s  refered  to  Zacks  and  Haber  [14].  Total  force  forecasting  procedures 
were  applied  by  the  Marine  Corps  following  the  technical  paper  of  Zacks 
[10].  The  method  was  also  applied  to  a problem  of  forecasting  the  reten- 
tion of  Navy  Pilots. 

4.  Strategies  of  Crossing  Mine  Fields. 

The  study  on  strategies  of  crossing  mine  fields  was  performed  in 
1963-1966  (see  Zacks  and  Goldfarb  [11]  and  Zacks  [8])  and  renewed  recently 
in  a new  study,  Zacks  [15].  The  problem  can  be  described  in  the  following 
terms. 

A number  of  clusters,  k,  containing  N^,...,N  mines,  are  distributed 
randomly  over  a (rectangular)  field  of  specified  dimensions  • 

The  centers  of  the  clusters  (C^.n^) ^k’\^  are  given  anci  ttie  sPeci“ 

fic  bivariate  distribution  functions,  according  to  which  the  mines  in  each 
cluster  are  distributed,  are  given  too.  The  problem  is  to  determine,  for 
any  specific  breaching  path,  the  survival  probabilities  of  each  one  of  n 
targets  (personnel,  vehicles  or  tanks)  crossing  in  a column  at  the  same  path; 
under  the  following  assumptions: 

(i)  the  encounters  of  targets  and  mines  are  conditionally 
independent  of  the  prior  events; 

(ii)  with  probability  a mine  is  detected  and  destroyed 

by  the  target; 

(iii)  an  activated  mine  destroys  the  target  only  with  proba- 
bility pk  ; 

(iv)  a target  passing  the  neighborhood  of  a mine  activates 
it  with  probability  p . 

The  parameters  Pdet»  Pk  and  Pact  are  sPeci^ic  to  the  type  of  target 
crossing  the  field  and  the  type  of  mine  used. 


f ~ - 

f ' — ■ ■■  -—■■■-■  jr— 7 
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The  survival  probabilities  of  the  targets  crossing  in  a column  de- 
pend in  addition  on  the  number  of  mines  in  the  path.  This  is  an  unknown 
variable,  whose  distribution  can  be  often  approximated  by  a Poisson  distri- 
bution, with  parameter,  X , which  depends  on  the  location  of  the  path 
relative  to  the  centers  of  the  k clusters;  on  the  number  of  mines 
N^,...,Njt  ; and  on  the  particular  distributions  of  the  clusters.  In  a 
recent  paper  Zacks  [15]  provided  formulae  for  the  determination  of  the 
parameter  X when  the  cluster  distributions  are  bivariate  normal,  and  an 
algorithm  for  the  recursive  (exact)  determination  of  the  survival  proba- 
bilities. In  the  previous  papers  of  Zacks  and  Goldfarb  [11]  and  of  Parsons 
[ 7 ] closed  form  formulae  were  provided  for  the  survival  probabilities  in 

the  special  case  of  no  possible  detection  (p,  = 0)  and  every  activation 

det 

destroys  the  target  (p^  = 1).  Closed  analytic  formulae  can  be  derived  also 
in  the  present  general  model.  They  are  however,  of  no  special  interest, 
because  the  computation  time  required  is  trivial.  The  procedure  developed 
by  Zacks  in  [ 15]  has  been  already  applied  by  TRASANA  at  the  White  Sands 
Missile  Base  for  the  solution  of  various  problems  in  the  evaluation  of  new 
weapon  systems. 

5.  Example . 

Consider  a squared  mine  field  of  dimensions  200m  x 200m.  Nine  clusters 

of  N * 50  mines  are  distributed  over  the  field  around  the  center  points 

with  coordinates  £ = 50,  100,  150  m and  * 50,  100,  150  m (i=l,2,3). 

Each  cluster  is  distributed  according  to  a bivariate  normal  distribution 

2 2 2 2 

around  its  center  point  with  variances  o = 81.25  m , o = 43.75  m and 

x y 

correlation  p » .5447.  Such  a distribution  is  obtained  by  a 0 = 30° 
rotation  of  a bivariate  normal  distribution  with  o * 10m  , a = 5m  and 


p * 0. 


-12- 


i 


Y’  Y X' 


Fig.  2.  Ellipsoids  of  concentration  of 
scatter  distributions. 


We  consider  straight  breaching  paths  of  width  lm.  We  compute  first  the 
probability  , that  a mine  belonging  to  cluster  (i,j)  will  fall  in 
the  breaching  path.  The  distribution  of  the  number  of  mines  in  the  path  is 


3 3 

approximated  by  a Poisson  distribution  with  mean  1*  « 50  S 2 i|i  . 

i-1  j-1  ij 

Following  the  algorithm  developed  in  [15]  we  can  determine  the  expected 
number  of  survivors  in  a column  of  n.  In  the  following  table  we  provide 
the  values  of  X*  and  expected  survivors  in  20  paths  located  B m from 
the  center  of  the  field. 


' 
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Table  4.  Expected  number  of  survivors  in  columns  of 


n-l°,  Pdet-.4, 


Pact=-95‘  V 


7>  Pdud”'2, 


n-50. 


B 

X* 

E{X  } 
n 

-95. 

0.000 

10.000 

-85. 

0.003 

9.999 

-75. 

0.114 

9.953 

-65. 

1.331 

9.452 

-55. 

4.552 

8.127 

-45. 

4.552 

8.127 

-35. 

1.334 

9.451 

-25. 

0.228 

9.906 

-15. 

1.334 

9.451 

- 5. 

4.552 

8.127 

5. 

4.552 

8.127 

15. 

1.334 

9.451 

25. 

0.228 

9.906 

35. 

1.334 

9.451 

45. 

4.552 

8.127 

55. 

4.552 

8.127 

65. 

1.331 

9.452 

75. 

0.114 

9.953 

85. 

0.003 

9.999 

95. 

0.000 

10.000 

Average 

9.259 

There  are  several  interesting  statistical  problems  that  can  be  tackled 
now.  One  is  the  problem  of  optimal  strategies  of  crossing  a field  when  there 
is  no  complete  information  on  the  scatter  distributions  of  the  mines  in  each 


cluster  and  on  the  number  of  mines  in  each  cluster.  If  this  Information  is 
available  we  can  easily  find  the  path  for  which  the  expected  number  of  sur- 
viving targets,  out  of  n crossing  in  a column,  is  maximal.  One  can  impose 
on  this  maximization  problem  also  additional  constraints  of  the  battlefield. 
On  the  other  hand,  when  the  information  needed  to  compute  the  distribution 
of  the  number  of  mines  in  each  path  is  incomplete  the  decision  problem  be- 
comes that  of  choosing  from  m alternative  paths  one  having  the  most 
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favorable  distribution  of  the  number  of  mines  in  the  path.  The  decision  of 
how  to  cross  the  n targets  is  performed  in  a sequential  manner.  According 
to  some  prior  consideration  a corssing  path  is  chosen  for  the  first  target. 
The  result  of  the  first  crossing  attempt  (success  or  failure  and  the  points 
along  the  path  at  which  mines  were  destroyed)  leads  to  a reevaluation  of  the 
strategy  and  a crossing  path  for  the  second  target  is  then  chosen,  and  so  on. 
The  objective  is  to  find  a decision  rule  for  maximizing  the  expected  number 
of  survivors.  Zacks  [8]  studied  the  problem  previously  for  the  case  of  two 
crossing  paths,  in  terms  of  the  famous  TWO  ARMED  BANDIT  PROBLEM.  It  is  an 
interesting  and  quite  difficult  problem  that  should  be  studied.  Algorithms 
of  the  Dynamic  Programming  type  for  the  optimal  determination  of  the  cross- 
ing strategy  can  be  developed. 
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